AO*- A 0A3  738  CORNELL  UNIV  ITHACA  N  Y  DEPT  OF  THEORETICAL  AND  APR— ETC  F/G  21/2 

MATHEMATICAL  THEORY  OF  LAMINAR  COMBUSTION.  V.  UNSTEADY  BURNING  — CTC(U) 
MAR  80  J  D  BUCKMASTER'  6  S  LUOFORD  DAAG29-79-C-0121 

UNCLASS  I F I  ED  TR-107 _ _ ARQ-158B2. 12-M _ NL 


MICROCOPY  RESOLUTION  TEST  CHART 

NAJJOMl  BURfAU  CM  STANDARDS  1 963-.’* 


MATHEMATICAL  ..THEORY  OF^AMINAR  .COMBUSTION. V* 
Unsteady  Burning  of  &  Linear  Condensate  . 


Technical  Report  No.  107 
J. D.  Buekmaster  &  G.S.S.  Ludford 
March  1980 


U.S.  Army  Research  Office 
Research  Triangle  Park,  NC  27709 

Contract  No.  DAAG29-79-C-0121 

Cornell  University 
Ithaca,  NY  14853 


Approved  for  public  release;  distribution  unlimited 


The  findings  of  this  report  are  not  to  be  construed 
as  an  official  Department  of  the  Army  position  unless 
so  designated  by  other  authorized  documents. 


Chapter  V 

Unsteady  Burning  of  a  Linear  Condensate 

1.  The  Problem 

Chapter  IV  was  concerned  with  the  steady  combustion  of  the  gases 
produced  by  vaporization  of  a  linear  condensate  at  its  surface,  whether 
•through  pyrolysis,  sublimation  or  evaporation.  The  results  were  character¬ 
ized  by  response  curves  of  burning  rate  versus  pressure  (represented  by 
the  Damkohler  number).  If  the  applied  pressure  varies  in  time,  then  so 
also  must  the  burning  rate;  it  is  the  nature  of  their  dependence  that  is 
examined  in  the  present  chapter. 

The  effect  of  variations  in  pressure  on  solid  pyrolysis  has  received 
a  considerable  amount  of  attention  because  of  its  relevance  to  the  stability 
of  solid-propellant  rocket  motors.  The  essential  idea  is  that  acoustic 
waves  bouncing  around  the  combustion  chamber  will  impinge  on  the  propellant 
surface  and  thereby  generate  fluctuations  in  the  burning  rate.  These 
fluctuations  will  affect  the  reflected  wave  which,  it  is  argued,  can  have 
a  larger  amplitude  than  the  incident  wave.  If  so,  the  transfer  of  energy 
(provided  it  is  greater  than  losses  through  dissipation  and  other  mechanisms) 
implies  instability. 

The  response  of  a  burning  condensate  (solid  or  liquid)  to  an  imping¬ 
ing  acoustic  wave  will  be  the  focus  of  our  discussion.  Mathematically  we 
must  deal  with  the  disturbance  of  a  steady  field  containing  large  gradients, 
so  that  a  frontal  attack  on  the  governing  equations  is  not  feasible.  Six 
regions  can  be  distinguished:  condensate ,. preheat  zone,  flame,  burnt  gas, 
entropy  zone  and  far  field;  so  that,  without  rational  approximation,  the 
discussion  soon  disappears  into  numerical  and/or  ad  hoc  analysis.  There 
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This  report  is  Chapter  V  of  the  twelve  in  a  forthcoming  research 
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are  three  parameters  other  thar  B~  1  whose  smallness  can  be  exploited, 
namely  the  Mach  number,  the  ratio  of  characteristic  time  in  the  gas  to  that 
in  the  condensate  and  the  ratio  of  densities.  Low  Mach  number  enables 
the  combustion  approximation  (with  its  spatial  constancy  of  pressure)  to 
be  used  in  the  combustion  regions,  namely  preheat  zone,  flame  and  burnt 
gas;  small  time  ratio  ensures  the  gas  phase  responds  much  more  rapidly 
than  the  condensate,  so  that  the  three  combustion  regions  are  quasi-steady 
if  the  acoustic  frequency  is  sufficiently  small  (a  restriction  of  no  practical 
significance) ;  and  small  density  ratio  implies  that  fluctuations  in  the 
location  of  the  surface  of  the  condensate  may  be  neglected.  The  activation- 
energy  asymptotics,  needed  for  an  analytical  description  of  the  combustion 
regions,  determine  how  small  the  acoustic  frequency  must  be:  the  resulting 
changes  in  flame  temperature  (due  entirely  to  fluctuations  in  the  con¬ 
densate)  must  be  0(0  ^) . 

Only  adiabatic  burning  will  be  considered;  no  work  on  heat  loss  has  so 
far  been  reported.  The  first  task  is  to  determine  the  response  of  the  burn¬ 
ing  rate  to  a  general  (small)  change  in  pressure  level  of  the  combustion 
regions,  without  regard  to  its  cause.  Sec.  2  determines  the  response  for 
solid  pyrolysis  and  Sec.  3  that  for  liquid  evaporation.  The  incident  acoustic 
wave  excites  such  a  response,  and  the  response  induces  a  reflected  wave, 
through  the  intermediary  of  an  entropy  zone  which  converts  isentropic  con¬ 
ditions  in  the  far  field  into  isothermal  ones  in  the  near  field  (Sec.  4). 
Conditions  for  the  reflected  wave  to  be  stronger  than  the  incident  wave  are 
developed  in  Sec.  5. 

Our  analysis  also  enables  stability  characteristics  to  be  discussed, 
at  least  on  the  time  scale  adopted,  by  holding  the  pressure  level  fixed 
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(i.e.  eliminating  the  incident  acoustic  wave).  The  stability  analysis  of 
Ch.  Ill  is  thereby  extended  to  what  may  be  called  anchored  flames.  The 
extension  continues  in  Sec.  6  to  a  flameholder,  where  the  combustion  field 
need  no  longer  be  taken  quasi-steady  but  nevertheless  still  varies  on  a 
long  time  scale.  Only  in  this  last  section,  where  the  combustion  field 
is  truly  unsteady,  does  the  Lewis  number  play  a  role. 

In  view  of  the  large  number  of  small  parameters  we  shall  eschew  writing 
formal  expansions.  Instead  the  sympol  <<  will  be  used  to  indicate  that 
only  the  leading  term  in  the  corresponding  parameter  is  being  considered. 

2.  Solid  Pyrolysis 

The  unsteadiness  will  be  treated  as  the  perturbation  of  a  steady 
state  with  Damkohler  number  Dq.  Then,  as  in  Sec.  III. 3,  the  burning  rate 
M0  in  the  steady  state  will  be  taken  as  the  mass  flux  on  which  units  are 
based;  according  to  the  result  (III. 16),  we  have 

Mo  "  T*exp(-e/2  Tj/0  (1) 


since  none  of  the  products  of  the  gaseous  reaction  is  produced  or  absorbed 
at  the  surface  of  the  condensate  (Js  ■  1) .  The  flame  temperature  T^  is 
to  be  calculated  from  the  formula  (IV. 3). 

The  smallness  of  the  Mach  number  ensures  that  the  combustion  approxima¬ 
tion  is  valid.  If,  in  addition,  the  fluctuations  in  the  applied  pressure 
are  not  too  rapid  then  the  combustion  is  quasi-steady  and  the  analysis  of 
Ch.  II  is  applicable.  More  precisely,  if  u  is  a  characteristic  frequency 

of  the  fluctuations  (e.g.  the  frequency  of  an  impinging  acoustic  wave)  then 

-1  2 
the  time  u  should  be  long  compared  to  the  response  time  pcX/cpMo  of 


the  gas  phase,  determined  by  diffusion.  In  3hort,  we  require 


a  restriction  of  no  practical  consequence  since  the  diffusion  time  is 
always  very  small.  The  same  formula  (1),  with  Mq  and  Dq  replaced  by 
M  and  D,  therefore  holds  for  the  burning  rate,  where  now  both  D  and 
are  functions  of  time,  the  former  prescribed  and  the  latter  to  be 
determined  from  an  analysis  of  heat  transfer  in  the  condensate  and  the 
nature  of  the  vaporization. 

This  conclusion,  which  is  not  restricted  to  solid  pyrolysis,  was 
reached  by  Denison  and  Baum  (1961)  using  ad  hoc  arguments  which  Williams 
(1973)  later  replaced  with  rational  asymptotic  analysis.  However,  none 
of  these  authors  noted  that  (in  general)  a  further  restriction  of  to  is 
required  if  the  Mach  number  is  to  remain  sufficiently  small  for  the 
combustion  approximation  to  apply.  In  our  formulation  the  restriction 
comes  from  allowing  the  flame  temperature  to  be  perturbed  only  by  0(0  *) 
amounts,  so  that  the  change  in  burning  rate  is  0(1).  More  precisely. 


M 


_l/2 
D  e 


♦  /  2 


(3) 


2 

where  Tjfr^/0  is  the  perturbation  (T^  now  being  the  unperturbed  flame 
temperature)  and  M,  D  are  measured  in  units  of  Mq,  Dq.  [The  formula 
follows  from  the  result  (III. 16).] 

Now  this  perturbation  is  due  to  fluctuations  of  temperature  in  the 
condensate,  and  it  may  be  determined  by  calculating  the  total  change  in 
enthalpy  they  produce  up  to  the  surface.  In  general,  slow  variations  on 

A  2  A 

the  time  scale  c  M  /p£  of  the  condensate  are  needed  if  the  result  is  to 
P  o 


-  5  - 


tt-a^/axe  .  (A) 

P  O 

Note  Chat,  since  the  condensate  is  much  denser  and  more  conductive  than 
the  gas,  i.e. 

Pc/p  <<  1  and  A/A  <<  1  ,  (5) 

the  change  in  enthalpy  in  the  gas  phase  is  negligible,  i.e.  in  response 
to  changing  conditions  the  gas  (unlike  the  condensate)  has  no  inertia. 

Thus  the  restriction  (4)  implies 

u  «  c  M^/p  xe  ,  (6) 

p  o  c 

whereas  rough  equality  would  have  to  hold  for  there  to  be  a  perturbation. 
[Note  that  the  condition  (2)  is  satisfied  a  fortiori.]  In  special  circum¬ 
stances  there  is  an  alternative  to  slow  variation,  to  which  we  shall  come 
later. 

To  calculate  the  change  in  enthalpy  we  consider  the  dimensionless 
temperature  equation 

3T/3t+M  3T/3x- 32T/3x2  -  0  (7) 

in  the  condensate.  Here  pA/c  M2  and  A/c  M  have  been  taken  as  units 

p  o  p  o 

of  time  and  length,  while  T  is  still  referred  to  Q/Cp.  Since  temporal 
variations  occur  on  the  long-time  scale 

T  -  t/0  ,  (8) 

the  term  3T/3t  ■  0  Ht/3t,  which  represents  local  storage  of  enthalpy, 
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ia  seen  to  be  a  perturbation.  If  it  is  neglected,  the  leading  approxima¬ 
tion 


T  »  T  +  (T  -T  )eix  (9) 

is  obtained,  where  M  varies  with  r.  The  analysis  so  far  is  valid 
whatever  kind  of  vaporization  occurs  at  the  surface  of  the  condensate; 
for  solid  pyrolysis,  T  is  connected  to  M  by  the  law  (IV. 5),  which 
now  reads 

M  -  (k/Mo)Texp(-0/Ts)  .  (10) 

From  the  leading  approximation  we  can  calculate  the  neglected  perturbation 
term,  which  represents  a  heat  source/sink  in  equation  (7),  and  hence  the 
perturbation  in  flame  temperature.  We  find 

°  3T  —9 

M*«  -  -/  fr  dx  -  -(p+q)M  ^dM/dx  ,  (11) 


where 


p  -  T?/(0+T  )  >0  and  q  -  T  -T  ,  (12) 

so  that,  from  the  result  (3  ),  the  equation  connecting  the  variations  in 
M  and  D  is 

* 

(p+q)dM/dx  -  M3iln(D/M2)  .  (13) 

The  steady  state  M  »  D  *  1  clearly  satisfies  the  last  equation;  and 
we  are  mainly  concerned  with  relating  the  small  disturbances  of  M  and  D 
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caused  by  an  acoustic  wave.  However  the  equation  may  also  be  used  to 
investigate  the  inherent  stability  of  the  steady  rate,  i.e.  the  way  M 
changes  when  D  is  held  fixed  at  1.  We  see  that  when  M  is  disturbed 
from  1,  either  up  or  down,  it  is  driven  further  away  whenever 

p  +  q  <  0,  i.e.  T  <  0T  /(0+T  )  .  (14) 

The  condition  is  difficult  to  interpret  because  Ts  itself  depends  on  M; 

but  certainly  there  is  linear  instability  if  it  is  satisfied  by  the  value 

of  T  corresponding  to  M  *  1  in  the  pyrolysis  law  (10).  This  instability 
s 

was  predicted  by  Denison  and  Baum  but  it  has  apparently  never  been  observed. 
Their  condition,  after  extraneous  terms  are  eliminated  by  letting 
0  -*■  «>,  is  actually 

T  <  0T  /(0-T  )  (15) 

s  s 

because  their  pyrolysis  law  lacks  the  factor  T  .  Figure  1  shows  the 

s 

stability  boundaries  for  a  factor  Tg  with  a  =  0,  1/2  and  1. 

For  a  general  pyrolysis  law  M(Tg)  the  right  side  of  the  inequality 

becomes  T  -  MdT  /dM,  and  it  is  conceivable  that  no  propellant  satisfies 
s  s 

such  a  condition.  The  procedure  for  checking  a  given  T  would  be  to 
compute  the  burning  rate  (1)  for  the  assumed  Dq  and  then  determine  Tg 
and  dTg/dM  for  M  =  1  from  the  pyrolysis  law. 

For  oscillations  about  the  steady  state  forced  by  an  acoustic  wave, 
we  write  * 

M  -  1  +  meia>T  ,  D  -  l  +  de1<UT  with  5  -  0u  (16) 


and  linearize.  Then 


m  *  cd  with  c  ^  2  +  itu(p+q)  (17) 

where  T  now  has  its  value  for  M  =  1.  This  relation  determines  the 
s 

reflected  wave  in  the  far  field,  as  we  shall  see  in  Sec.  4.  [The 
.instability  result  above  is  obtained  by  requiring  c  ^  to  vanish.] 

An  alternative  way  for  changes  in  flame  temperature  to  be  0(0  is 
for  the  temperature  in  the  solid  to  fluctuate  by  that  amount  everywhere. 
The  fluctuations  can  then  be  on  the  scale  of  t,  rather  than  t,  so  that 
we  are  faced  with  solving  the  full  heat  equation  (7)  with  M  an  arbitrary 
function  of  t.  A  general  analysis  would  therefore  be  quite  complicated; 
accordingly  just  the  case  of  interest  will  be  considered,  namely  small 
departures  from  the  steady  state. 

When 


M  -  l  +  meluit  and  D  »  l  +  de1^  with  m,  d  «  1  (18) 

the  temperature  in  the  solid  has  a  steady  component,  given  by  the  distribu¬ 
tion  (9)  with  M  =  1  and  T  the  corresponding  surface  temperature,  and 

s 

a  fluctuating  component 

m[(-iq/aj)exp(x+iut)  +  (p+iq/w)exp(icx+iu>t) ]  with  <  *  (1  +  /l+4ioj)/2 

(19) 

The  two  contributions  to  this  result  are  due  to  fluctuations  in  M  within 
the  solid  and  at  its  surface  respectively.  We  now  find 

0  *  -/  -P  dx  ■  (1-k)  (p+q/K)meltllt  , 


(20) 


in  place  of  the  result  (11),  so  that 

m  *  cd  with  c  *  =  2  -  9(l-tc)  (p+q/nc)  ,  (21) 

which  reduces  to  the  previous  response  (17)  for  small  but  is  quite 
different  otherwise,  as  we  shall  see  in  Sec.  4. 

The  inherent  stability  of  the  steady  state  is  determined  by  setting 
c  ^  *  0.  Then 

Pk2+  (Q-P+2)k-Q  =  0  with  P  =  9p,  Q  =  9q  ;  (22) 

for  instability  the  parameter  values  P,  Q  must  make  the  real  part  of 
iu>  *  k(k-1)  positive  for  a  k  whose  real  part  is  positive.  Fig.  2  shows 
the  region  of  instability  in  the  P,  0-plane;  its  boundary  asymptotes 
P  +  Q  =  -6,  in  accord  with  the  previous  result  (14)  for  P,  Q  large. 

For  consistency,  we  must  have 

p,q  -  0(9-1)  ,  (23) 

which  certainly  requires  both  Ta  and  T  ^  to  be  of  the  same  order. 

Such  small  temperatures  do  not  mean  that  the  condensate  is  cold,  but  only 
that  the  corresponding  enthalpies  are  small  compared  to  the  heat  of 
combustion,  which  is  used  as  unit.  The  requirement  (IV. 6)  for  the  steady 
state  does  however  mean  that  Tw,  as  calculated  from  the  formula  (IV. 3), 
must  be  no  greater  than  1  (i.e.  the  pyrolysis  must  be  endothermic). 

3.  Liquid  Evaporation 

Liquid  propellants  have  also  been  used  in  rocket  motors,  though  not 
in  the  manner  envisaged  here.  A  better  motivation,  if  one  is  needed,  is 
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a  solid  propellant  which  first  liquefies  at  its  surface  and  then  vaporizes. 
If  attention  is  again  focused  on  adiabatic  combustion,  the  analysis  follows 
that  in  the  last  section  until  the  pyrolysis  law  is  reached,  when  the 
Clausius-Clapeyron  relation  (IV. 21)  is  used  instead,  i.e.  we  write 

D  -  kTe(T  -T  )-1exp (-0/T  )  .  (24) 

S  “  s  s 

Here  we  have  taken  the  Damkohler  number  proportional  to  p^  in  accordance 
with  the  formula  (1.59)  for  a  first-order  reaction  (otherwise  a  power  of 
D  is  required);  and  various  parameters  have  been  absorbed  into  k. 

Since  Tg  is  now  a  function  of  D,  the  latter  enters  into  the  temper¬ 
ature  distribution  (9)  in  the  liquid,  so  that  the  formula  (11)  for  the 
change  in  flame  temperature  becomes 

M<|>  -  -(qM‘2dM/dT  +  pM“1D~1dD/dT)  (25) 

OO 

where  now 

P  “  T2(T  -T  )/[T2+(BT  +9)(T  -T  )]  >  0  .  (26) 

8  w  S  S  S  00  S 

The  equation  connecting  the  variations  in  M  and  D  is  therefore 
qM-1dM/dT+pD-1dD/dT  -  M2in(D/M2)  . 

In  particular,  there  is  inherent  instability  (D  =  1)  when  q  is  negative, 
i.e.  for 


(27) 


Note  that  this  instability  (which  has  not  been  reported  before)  is 
independent  of  the  precise  form  of  the  vaporization  law  at  the  surface: 
it  occurs  whenever  the  surface  temperature  is  determined  by  the  surface 


pressure. 

For  general  small  departures  (18)  from  the  steady  state  we  find 

m  =  cd  with  c  *  (l-iip)/(2+iiq)  .  (28 

This  result  will  be  needed  when  we  come  to  acoustic  response  in  Sec.  5. 
When  p  and  q  are  0(8  ^)  higher  frequencies  w  can  be  considered, 
as  for  solid  pyrolysis.  The  fluctuating  component  (19)  in  the  condensate 
now  becomes 

(— iqm/ w)  exp (x+iut)  +  (pd+iqm/o))exp(KX+iut)  , 

where  p  has  its  new  definition  (26),  so  that 

m  »  cd  with  c  *  [1  +  0 (1— ic)p] /  [2  -  6(l-ic)q/ic]  .  (29^ 

When  D  is  held  fixed  at  1,  i.e.  d  “  0,  the  denominator  of  c  must 
vanish,  so  that  k  -  Q/(2+Q).  Instability  requires  k  >  1,  i.e. 

Q  <  -2  with  Q  -  9q  ,  (30 

a  condition  in  accord  with  the  earlier  result  (27)  for  Q  large. 

A.  Response  to  an  Impinging  Acoustic  Wave 

Analysis  of  the  combustion  zone  provides,  amongst  other  things,  a 
complete  description  of  the  state  of  the  burnt  gas  behind  the  flame  sheet 
in  terms  of  the  time-dependent,  but  spatially  uniform,  pressure  there. 
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Conditions  are  essentially  isothermal  with  the  velocity  determined 

by  the  instantaneous  value  of  M. 

The  varying  pressure  is  generated  by  some  source  far  from  the  flame; 
we  shall  confine  ourselves  to  acoustic  waves  travelling  normal  to  the 
surface  of  the  condensate,  which  immediately  raises  a  problem.  The  far 
field  is  isentropic  when  the  near  field  is  isothermal,  so  that  the  two 
cannot  be  matched.  In  other  words,  it  is  not  simply  a  matter  of  evaluating 
the  pressure  in  the  far  field  as  the  near  field,  represented  by  a  single 
point,  is  approached.  An  "entropy"  layer  is  needed  to  transform  con¬ 
ditions  from  isentropic  to  isothermal  in  a  distance  small  on  the  scale 
of  the  far  field  but  large  on  that  of  the  near  field. 

In  the  chemistry-free  region  beyond  the  flame  sheet,  the  equations 
with  which  we  have  to  deal  are  those  of  a  compressible,  heat-conducting, 
viscous  fluid,  i.e.  the  dimensional  form  of  equations  (1.53,  55,  56)  with 
pressure  and  viscous  terms  restored  in  the  energy  balance,  with  £2  =  0 
and  with  the  perfect-gas  law  reinstated.  Linearized  about  the  steady 
state  pa,  V",  pc,  Tw  they  read 

3p/3t  +  voo3p/3x  +  po#3v/3x  =  0  ,  (31) 

p^Ov/St  +  v^Sv/Sx)  -  -3p/3x+  (4</3)32v/3x2  ,  (32) 

p  c  (3T/3t  +  v  3T/3x)  -  X32T/3x2  -  3p/3t  +  v3p/3x  .  (33) 

00  p  00  00 

It  is  natural  to  take  p  .  a_  -  /yRT  /m,  p  -  p  aVy  and  T  as  units 

OOOD  00  C  00  00  00 

of  density,  velocity,  pressure  and  temperature  in  such  a  combustion-free 
region.  Length  will  be  referred  to  a  distance  Z  to  be  specified  later 
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and  time  to  A/a  .  Then,  for  disturbances  proportional  to  e*Ut,  the  dimension- 
00 

less  equations  are 

-1  2  2 

imp  +  m 3p/3x  +  3v/3x  -  0  ,  iwv  +  m  3v/3x  *>  -y  3p/3x+KD3  v/3x  , 

00  Jt, 

(34) 

iuT  +  m^ST/Sx- X£32T/3x2  -  (y-1)  (iup +  mo3p/3x)/y  ,  p  -  T  +  p 

(35) 

where 

-  bK/3poaa{ol  and  X£  =  X/CpP^A  .  (36) 

Solution  of  these  equations  determines  the  dimensional  density,  velocity, 
pressure  and  temperature  as 

P00(l+peiut)  ,  +  a^ve*1^  ,  pc(l+peiwt)  and  Ta>(l+Telu>t)  (37) 

respectively. 

These  equations,  being  linear  with  constant  coefficients,  can  be 
solved  exactly  and  the  solution  used  to  trace  the  transition  from  isothermal 
to  isentropic  conditions.  However,  to  see  how  the  burning-rate  fluctuations 
are  caused  by  the  incident  acoustic  wave  and  then  how  they  determine  the 
reflected  wave,  across  the  entropy  layer,  it  is  easier  to  use  expansions 
appropriate  to  the  various  regions. 

To  describe  the  acoustic  waves,  we  choose  for  l  the  wavelength 

*a  “  a»/u  >:>  *c/m-  >:>  *c  <38) 

where  *c  -  X/CpP^  is.  to  within  a  factor  pj pB,  the  unit  of  length 
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used  in  Che  combustion  zone.  We  have  used  the  restriction  (2) ,  which  now 
reads 


id  «  v ,  (39) 

oo  c 

and  the  fact  that  the  Mach  number  m  is  small  to  obtain  these  extreme 

00 

inequalities.  It  follows  that  the  dimensionless  to  is  1  but  that 
and  A^  are  small  so  that 

ix  -ix 

p  =  Aea  +  Be  3  ,  T  =  (y-l)p/Y  »  (40) 

ix  -ix 

P  ■  p/Y  >  v  -  (-Ae  a  +  Be  a)/y  (41) 

where  x  is  distance  measured  in  units  of  £  .  Here  A  and  B(m  ) 

3i  d  oo 

are  the  amplitudes  of  the  pressure  in  the  incident  and  reflected  waves, 
respectively,  which  would  be  dissipated  over  a'  distance 

L  ■  m  £  f\  »  £  /m  »  £  (k  is  A  for  £  *  £  )  .  (42) 

a  ®aa  a®  a  a  »  a 

For  the  theory  to  be  applicable  Lfl  must  be  large  compared  to  an  overall 
dimension  but,  since  the  latter  is  normally  comparable  to  £&,  that 
provides  no  restriction  in  practice.  The  goal  is  to  determine  B  when 
A  is  given,  so  as  to  calculate  the  reflection  coefficient  [b/a|. 

On  the  acoustic  scale  the  whole  combustion  field  is  represented  by 
a  single  point,  which  we  may  take  to  be  xft  *  0.  The  temperature  disturbance 
there  is  zero,  which  clearly  contradicts  the  isentropic  requirement  (40b) 
of  the  acoustic  waves  when  there  is  a  pressure  disturbance.  A  layer  is 
needed  to  change  the  entropy  from  its  constant  value  in  the  acoustic  field 
at  0.  Now,  in  the  absence  of  dissipation,  entropy  is  convected 


with  the  fluid;  to  describe  these  waves,  we  choose  for  SL  the  wavelength 


l 

e 


v  /ii>  =  mi  «  i, 

oo  od  a  a 


(A3) 


The  dimensionless  frequency  equals  m^  and  hence  is  small.  We  also 
find  that  and  are  still  small,  their  dissipative  effect  being 

felt  over  a  distance 

\  “  m.VXe  >>  *e  (Ae  ls  Xi  for  1  =  V  *  (44) 


If  is  to  be  small  compared  to  we  must  have 


u  »  m  v  /£ 

00  CO 


c 


f 


(A5) 


which  now  places  a  lower  bound  on  the  frequency  for  the  theory  to  apply. 
The  restriction  is  of  no  consequence  in  practice  but,  if  it  is  not 
introduced,  entropy  effects  extend  into  the  far  field  and  our  simple 
picture  is  destroyed.  The  results  must  then  be  obtained  from  the  full 
solution  of  equations  (34,35). 

We  now  come  to  the  structure  of  the  entropy  layer,  which  the  momentum 
equation  shows  to  be  isobarlc,  as  might  have  been  expected.  The  complete 
solution  may  therefore  be  written 

T  =  (Y-l)p[l-exp(-ixe-Xe)]/y  ,  p  -  p[l+(y-l)exp(-ixe-Xe)]/y  , 

v  »  moe[-ipxe/Y+  (c-l)p]  (46) 

to  leading  order,  where  x.  and  X  are  distances  measured  in  units  of 

e  e 

Ae  and  Lg,  respectively.  Here  the  constants  of  integration  have  been 


determined  so  that 
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T  ■  0  and  p  +  v/m<o  »  cp  (47) 

at  x  ■  X  =  0,  as  required  by  the  analysis  of  the  combustion  regions:  to 
6  € 

0(1)  accuracy  in  9  the  temperature  is  unaffected  and  m  *  cd  (in  the 
notation  of  Secs.  2  and  3).  Note  that  v  as  xg  -*■  «®,  a  fact  that 

makes  a  careful  analysis  of  the  entropy  layer  necessary,  as  we  shall  see 
shortly. 

As  Xg  -*■  00  the  solution  (46)  becomes  isentropic,  i.e.  the  relations 
(40b,  41a)  are  satisfied.  A  formal  matching  is  obtained  by  taking  m^ 
as  the  vanishingly  small  parameter,  which  is  suggested  by  the  transforma¬ 
tion  x  ■  m  x  .  When  applied  to  v  the  matching  shows  that 
d  ®  6 

B  ■  A(l+m  a)  with  a  ■  2y(c-l)  (48) 

QO 

correct  to  first  order  in  m^.  The  reflected  wave  is  thereby  determined 

in  terms  of  the  incident  wave.  tUse  of  m^  as  small  parameter  makes 

3  2 

the  inequalities  (39)  and  (45)  become  0(11^)  <<  to  <<  000.] 

Usually  such  reflection  problems  are  discussed  in  terms  of  the 
acoustic  admittance.  In  the  present  context  Williams  (1973)  equates  the 
value  of  v/p  as  x  -►  0  to  its  value  in  the  combustion  region,  which 
happens  to  give  the  correct  result  (48).  It  would  be  more  logical  to 
equate  it  to  the  value  at  the  edge  of  the  entropy  layer,  but  that  is 
impossible  since  v  ->■  »  as  x  -»■  «.  [The  difficulty  is  not  a  creature 
of  our  analysis:  the  continuity  equation  (31)  requires  v  to  become 
unbounded  whenever  p  does  not  vanish  at  the  edge  of  the  layer.]  Only 
proper  matching  of  the  entropy  layer  resolves  the  delemna. 


■  k.*. 


r 
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5.  Amplification 

Since  A  is  the  amplitude  of  the  incident  wave  and  B  is  that  of 
the  reflected  wave,  there  is  amplification  for  small  values  of  db  if 
and  only  if 

Re(a)  >  0  i.e.  Re(c)  >1  .  (49) 

We  shall  now  examine  the  two  cases  of  solid  pyrolysis  in  Sec.  2  and  the 
two  cases  of  liquid  evaporation  in  Sec.  3  for  this  instability. 

The  response  (17)  for  solid  pyrolysis  shows  what  Re(c)  is  never 
greater  than  1/2,  its  value  for  u  =  0,  so  that  there  is  no  instability. 
The  conclusion  is  changed  by  giving  D  a  pressure  dependence  pn  with 
n  >  2  since  then  c  is  replaced  by  nc,  whose  maximum  real  part  exceeds 
1;  but  in  practice  n  is  never  that  large. 

On  the  other  hand  the  response  (21)  shows'  that 

Re(c)  *  1/2+ [2P  +  4- (QfP+2)2]u)2/8  +  0(u4)  (50) 

for  small  to,  which  increases  initially  with  u  for  a  band  of  values 
around 


Q  -  -(P+2)  ,  (51) 

suggesting  that  the  inequality  (49)  may  be  satisfied  for  some  w  when 
P  and  Q  are  chosen  appropriately.  This  is  confirmed  by  graphs  of  Re(c) 
versus  u  for  several  values  of  P  in  Fig.  3,  where  the  value  (51)  has 
been  taken  for  Q  to  ensure  the  fastest  rate  of  initial  increase. 

Turning  now  to  the  evaporating  liquid,  we  find  that  the  response 


(28)  gives 


For  Instability  q  must  be  negative  but  greater  than  -p;  then  all 
frequencies 

u>  >  /-2/q(p+q)  (53) 

are  amplified.  On  the  other  hand,  Re(c)  for  the  response  (29)  behaves 
like  -P/w//2(2+Q)  as  ai  -*•«»,  so  that  1  is  certainly  exceeded  for  Q  <  -2 
whatever  the  value  of  P.  Figure  A  shows  graphs  for  P  •  1  and  various 
values  about  Q  ■  -2. 

The  absence  of  amplification  at  all  frequencies  to  for  a  pair  of 
values  p,  q  does  not  necessarily  mean  that  stable  burning  will  take 
place.  There  may  be  inherent  instability.  For  solid  pyrolysis,  para¬ 
meter  values  for  which  p  +  q  is  negative  are  ruled  out  even  though  there 
is  no  amplification  of  acoustic  waves.  For  liquid  evaporation  there  is 
stable  burning  only  when  q  is  positive.  Similar  remarks  apply  to  the 
cases  governed  by  P  and  Q. 


6.  Stability  of  Anchored  Flames 

The  flames  considered  in  this  chapter  are  anchored  to  the  condensate 
in  the  sense  that  they  are  controlled  by  the  manner  in  which  the  reactant 
is  liberated  at  the  surface  (which  in  turn  is  affected  by  heat  from  the 
flame).  The  stability  of  certain  anchored  flames  has  therefore  been 
covered  already.  Here  we  shall  consider  a  stronger  control  in  which  there 
is  no  influence  of  the  flame  on  the  source,  i.e.  the  reactant  is  supplied 
at  constant  temperature  Tg,  any  heat  conducted  back  from  the  flame  being 
Instantly  removed.  The  flamholder  of  Sec.  II. A  is  an  example,  if  there 
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is  temperature  control  at  the  porous  plug  (by  means  of  heating/cooling 
coils)  and  the  flow  rate  M  is  maintained  constant.  While  stability  is 
the  main  interest,  we  shall  also  consider  the  effect  of  an  impinging 
acoustic  wave,  i.e.  oscillations  in  D. 
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